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Abstract. Recent measurements of elliptic flow (^2) and the nuclear modification 
factor {Rcp) of strange mesons and baryons in the intermediate px domain in Au+Au 
collisions demonstrate a scaling with the number-of-constituent-quarks. This suggests 
hadron production via quark coalescence from a thermalized parton system. Measuring 
the elliptic flow of charmed hadrons, which are believed to originate rather from 
fragmentation than from coalescence processes, might therefore change our view of 
hadron production in heavy ion collisions. 

While direct V2 measurements of charmed hadrons are currently not available, 
single electron V2 at sufficiently high transverse momenta can serve as a substitute. 
At transverse momenta above 2 GeV/c, the production of single electrons from non- 
photonic sources is expected to be dominated by the decay of charmed hadrons. 
Simulations show a strong correlation between the flow of the charmed hadrons and 
the flow of their decay electrons for pr > 2 GeV/c. 

We will present preliminary STAR results from our single electron V2 measurements 
from Au+Au collisions at RHIC energies. 
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1. Introduction 

Elliptic flow of hadrons has long been suggested as a sensitive probe of the pressure built 
up in the early stages of heavy ion collisions . It has furthermore been proposed as a 
signal of a phase transition from Quark-Gluon-Plasma to hadronic matter pi . At RHIC 
energies, for the first time in relativistic heavy ion collisions, elliptic flow of unflavored 
and strange hadronsj was found to be reasonably well described by hydro-dynamic 
calculations jSj up to transverse momenta of at least pt = 1-5 GeV/c as shown in Figure 
IH4j for the 80% most central events from Au + Au collisions at ^/snn = 200 GeV. 

At transverse momenta 
above pt = "2 GeV/c how- 
ever, the measured V2 val- 
ues seem to saturate and 
therefore are not describe by 
the hydro-dynamic calcula- 
tions anymore. More inter- 
esting than this deviation be- 
tween data and calculation is 
an apparent splitting between 
mesons and baryons, indepen- 
dent of their masses. While 
mass ordering is an essen- 
tial feature of hydro-dynamics 
and has been shown to lead 
to lower V2 values for heavier 
particles in the low transverse 
momentum region j^, a mass 
independent splitting between 
mesons and baryons has also 
been observed in the nuclear 
modification factor Raa [6 . Voloshin first suggested a scaling of V2 and Raa with the 
number of constituent quarks for all particle species 0, which has been interpreted 
as a signal for a partonic phase in the early stages of the collision. This interpreta- 
tion is supported by quark coalescence and quark recombination models, which are 
able to reproduce the measured f 2 and Raa values at intermediate transverse momenta 
{Pt ~ 2 — 6 Gev/c) when assuming a system of thermalized quarks 0. Please note, at 
even higher transverse momenta {pr > 6 GeV/c) , hadron production is expected to be 
dominated by parton fragmentation, hence quark coalesence models should break down 
in the high-p^ region ^9]. 

I Throughout this manuscript we wiU caU the u and d quarks the unflavored quarks. Hadrons built up 
exclusively from unflavored constituent quarks will be called unflavored hadrons. In this notation, the 
0-meson will be flavored, u, d and s quarks will be called light quarks. 
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Figure 1. Elliptic flow amplitude V2 for various 
strange hadrons as a function of transverse 
momentum pT- 
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In the regime where quark coalescence is vahd, the elhptic flow for all particle 
species can be expressed as 



V2h 



/I \ 

i=0 ^ ^ 



(1) 



Here, n is the number of constituent quarks and v^h {v^q^ is the hadron (quark) elliptic 
flow amplitude. 

Figure El |4j shows the scaled flow amplitudes W2/n as a function oi pt/^ for three 
meson and three baryon species, again for the 80% most central events from Au + Au 
collisions at y^J/vw = 200 GeV. Ignoring the low momentum pions, which suffer from 
strong resonance feed-down, the number-of-constituent-quark-scaling seems to nicely 
hold for hadrons built up from u, d and s quarks. Note, there is no evident ordering 
between hadrons of different strangeness content, suggesting that the s quark's flow is 
similar to that of the m, d quarks jH] • 
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Figure 2. Scaled elliptic flow amplitudes V2/n as a function oipr/n for mesons and 
baryons. 

With the arrival of the first charm measurements from the RHIC detectors, it is now 
of great interest to study the elliptic flow of charmed hadrons. Within quark coalescence 
models, and V2 of d and s quarks being determined from Figures Hand El elliptic flow 
measurements from D-mesons present a direct measure of the c quark elliptic flow (f2c)- 
In contrast to the m, d and s quarks which are predominantly produced in soft collisions, 
c quarks are expected to be produced via hard processes. At the time of their creation, 
c quarks should therefore be unaware of the collision geometry, i.e. carry no azimuthal 
anisotropics. Due to the large mass of the c quark, Dong et al [TT] argue, that large V2c 
values can only be achieved if re-scattering in a partonic phase is strong enough to also 
thermalized the light quarks. 

The measurement of D-meson elliptic flow might therefore not only challenge our 
current understanding of charmed hadron production, but also strengthen the case of 
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a system with partonic degrees of freedom and is a sensitive probe of thermalization 
therein. 

2. Concepts 

Although we have reported on measurements of open charm directly via invariant 
mass reconstruction of D-Mesons from hadronic decays, i.e. D ^ K -\- tt, in d+Au 
collisions ^U], low signal statistic and the large combinatoric background beneath the 
signal do not allow to study azimuthal anisotropies. However, at higher transverse 
momenta semileptonic decays {D e + X) are expected to dominate the single electron 
spectrum and for momenta higher than p > 2 GeV/c, the decay electron's azimuthal 
angle reflects the D-meson's azimuth angle very accurately as shown in the left panel of 
Figure El Plotted is the cosine of the difference of the D-meson's and its decay electron's 
azimuthal emission angles {(p Electron — (po) is versus the electron momentum. The dashed 
box reflects the momentum region used for this manuscript. Simulations done for this 
analysis and in [TTl |T31 EI suggests that single electron f 2 is a good measure of open 
charm flow as seen in the right panel of Figure El (taken from 12 ). Here, the solid 
lines depict the D-meson flow from quark coalescence under the assumptions that the 
c quark a) does not flow (t>2c = 0, lower line), and b) flows just as the light quarks 
{v2c = V2q-iight, upper line). The dashed lines and the symbols represent the flow of 
single electrons from D-meson decays under these assumptions. 
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Figure 3. (Color online) Left: Correlation between the difference of the D-meson's 
and its decay electrons's azimuth emission angles {4> Electron — 4>d) and the electron's 
momentum. Right: Elliptic flow of D-mesons (solid lines) and of single electrons from 
their decays (dashed hues, symbols). Figure taken from [T^. 



Note, for both assumptions, only D-meson production from coalescence is taken 
into account. Possible anisotropies from jet fragmentation are omitted. 
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3. Analysis 

For the analysis reported on here, we studied data from Au+Au collisions at y/s^N = 
200 GeV recorded by STAR^Hl during Brookhaven's RHIC beam-time in the year 2001. 
Selecting the 80% most central events from a minimum bias triggered data set and after 
applying quality cuts (e.g. requiring the primary vertex to be located within ±35 cm 
from the center of the STAR Time Protection Chamber (TPC) ^3] along the beam- 
axis), we used about 2M events. 

Similar as described in J7], the event-by-event reaction-plane was determined by 
calculating the Q- vector from charged particles within the phase space region 0.1 < 
p/ (GeV/c) < 2 and \ri\ < 0.8. For the electron candidate tracks to be correlated with the 
reaction-plane the momentum range was restricted to 2 < p/(GeV/c) < 3 and a cut on 
the track's specific ionization of the TPC gas was placed at 3.6 < dE / dx / {keV / cm) < 5. 
The choice of these selection criteria is due to the vr, K, p, and deuteron dE/dx bands 
crossing over the electron band at lower momenta and the limited separation between 
the bands above p = 2 GeV/c as seen in the left panel of Figure HI 




Figure 4. (Color online) Left: The specific ionisation of the TPC gas as a function 
of logio{p/ [GeV I c)). Right: Distribution of the specific ionisation of the TPC gas 
in the momentum region 2 < p/{GeV/c) < 3. The three gaussian curves show the 
contribution of tt, K, p and d centered at about 2.6 keV/cm, the contribution of 
electrons at about 3.5 keV/cm, and backgroimd from merged tracks at about 6 keV/cm. 

The right panel shows the dE/dx distribution in the momentum range 2 < 
p/(GeV/c) < 3. The three gaussian curves represent: a) the combined contributions 
from TT, K, p, and deuteron tracks centered at about 2.8 keV/cm, b) the contribution 
from electrons centered at ~ 3.5 keV/cm, and c) background believed from merged 
tracks§ at 6 keV/cm. The shaded areas indicate the dE/dx region accepted, leading 
to a non-electron contamination of about 10%, while rejecting about 80% of the total 

§ two tracks close together in the TPC which are falsely identified as one single track 
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electron yield. 

Before these electron candidates are correlated with the reaction-plane, we perform 
an invariant mass calculation for each candidate with all other tracks within the same 
event. The so obtained invariant mass distribution (see Figure Ej) features three distinct 
regions: a) a sharp peak from 7 — > e+ + e~ conversion in the detector material, b) 
a smeared out peak from the tt*' — > e"*" + e~ + 7 Dalitz decay where the 7 is not 
reconstructed, and c) the combinatoric background of random e+, e~ pairs. All electrons 
which form a pair with an invariant mass of less that 100 MeV/c^ are removed from the 
candidate sample. 

With detailed simulations of 
the STAR detector and using the 
TT^ distribution from ^H], we find 
that this method eliminates about 
50% of all electrons originating 
from tt" decays and 7-conversions. 
In the momentum interval 2 < 
p/(GeV/c) < 3, we estimate our 
electron sample to constitute to 
63% from D-meson decays and to 
37% from electrons from the re- 
maining photonic sources (7 con- 
version and 7r° Dalitz decays). Elec- 
tron backgrounds originating from 
other sources, e.g. the rj Dalitz de- 
cay, are expected to be relatively 
small and are neglected for this 
analysis. 

The azimuth angles from the 
remaining e^ are now correlated 
with the reaction-pane angles = 
atan2 {Qy, Qx) and the cos {2[(f) — ipRp]) distribution (see Figure IH)) recorded. The V2 
value can then be calculated V2 = {cos (2[(l) — ip^p])) /'ip^p . Here, ipRp ~ 0.7 is the 
reaction-plane resolution which is determined via the sub-event method as described in 

m- 

4. Results 

In order to correct for the 37% background from photonic sources, we again use the 
7r° spectrum from JH] onto which we impose elliptic flow according to Figure |21 [v^""^ = 
17%). Passing the resulting azimuthal anisotrop vr^-distribution through our detector 
simulation and analysis code, we obtain a cos (2[0 — ipRp]) distribution as shown by the 
squares in Figure El Subtracting this distribution from the uncorrected cos (2[0 — ipRp]) 
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Figure 5. Electron removal via e+e 
invariant mass method. 
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Figure 6. (Color online) cos (2[0 — ipRp]) distribution for uncorrected (solid bullets), 
corrected (open bullets) and background (squares) 



distribution (solid bullets), we obtain the cos {2[(j) — ipfip]) distribution for from 
D-meson decays (open bullets). Again, the elliptic flow amplitude is calculated as 
V2 = {cos (2[0 — ipRp])) /ip^p and evaluates to f2e± ~ 0.15±0.02(stat) in the momentum 
interval 2 < p/{GeV/c) < 3. Our preliminary estimate of systematic uncertainties 
originates from uncertainties in the V2n» and the exact amount of photonic e^ background 
and adds up to ~ 25% (relative). Other sources of systematic uncertainties (e.g. a slight 
systematic bias against central events when removing 7-conversion electrons) are still 
under investigation. 

Figure [3 shows our preliminary single electron elliptic flow results as a function 
of transverse momentum (open bullets). Again, the vertical bars represent only the 
statistical uncertainties and our preliminary systematic uncertainty estimate is about 
25%. The flgure also features data from a similar analysis presented by the Phenix 
experiment |E] (full bullets). Neglecting the difference in the centrality of the event 
samples used (Phenix: min. bias, this analysis 0-80% central), our data agree well with 
the Phenix results and extend them smoothly to higher momenta. 

The data points seem to favor the prediction from the quark coalescence model 
|121 IT^ under the assumptions of a partonic stage with thermalized and strongly 
flowing c quarks (solid line). However, with the present understanding of systematic 
uncertainties, the non-flow assumption can not be ruled out completely. 
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Figure 7. (Color online) V2 from non-photonic e± vs transverse momentum in Au+Au 
collisions at \/Nnn ~ 200 GeV. The open bullets represent preliminary data from 
this analysis, the full bullets represent results from 14 . The solid(dashed) line show 
the expectations from a quark coalescence model with thermalized and flowing(non- 
interaction and azimuthal isotrop) c quarks. 
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